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Abstract — Folding a sequence S into a multidimensional box 
is a method that is used to construct multidimensional codes. 
The well known operation of folding is generalized in a way 
that the sequence S can be folded into various shapes. The new 

definition of folding is based on lattice tiling and a direction in 

_i 

the D-dimensional grid. There are potentially 2 different 
folding operations. Necessary and sufficient conditions that a 
lattice combined with a direction define a folding are given. The 
immediate and most impressive application is some new lower 
bounds on the number of dots in two-dimensional synchroniza- 
tion patterns. This can be also generalized for multidimensional 
synchronization patterns. We show how folding can be used 
to construct multidimensional error-correcting codes and to 
generate multidimensional pseudo-random arrays. 

I. Introduction 

Multidimensional coding in general and two-dimensional 
coding in particular is a subject which attract lot of atten- 
tion in the last three decades. It includes error-correcting 
codes [1], [2], [3], [4], synchronization patterns [5], [6], 
[7], [8], [9], perfect maps and pseudo-random arrays [10], 
[11], and other topics as well. But, although the related 
theory of the one-dimensional case is well developed, the 
theory for the multidimensional case is developed rather 
slowly. This is due that the fact the most of the one- 
dimensional techniques are not generalized easily to higher 
dimensions. Hence, specific techniques have to be developed 
for multidimensional coding. One technique that was used 
for multidimensional coding is folding [8], [9], [10]. It was 
used to form a two-dimensional array, in the shape of a 
rectangle, from a one-dimensional sequence. It is generalized 
for multidimensional arrays to form a multidimensional box. 

In this paper we generalize the definition of folding. It is 
generalized in a way that we will be able to apply it to a 
multidimensional shape S which does not necessarily has a 
shape of a box. We present applications of the new definition 
to form multidimensional synchronization patterns, error- 
correcting codes, and pseudo-random arrays, for many types 
of multidimensional shapes. 

The rest of this paper is organized as follows. In SectionHJ 
we define the basic concepts of folding and lattice tiling. 
Tilings and lattices are basic combinatorial and algebraic 
structures. We will consider only integer lattice tilings. We 
will summarize the important properties of lattices and 
lattice tilings. In Section |lll] we will present the general- 
ization of folding into multidimensional shapes. We will 
show that all previous known foldings are special cases 
of the new definition. The new definition involves a lattice 
tiling and a direction. We will prove necessary and sufficient 



conditions that a lattice with a direction define a folding. In 
Section |IV] we give a short summary on synchronization 
patterns and present basic theorems concerning the bounds 
on the number of elements in such patterns. In Section [V] 
we apply the results of the previous sections to obtain new 
type of synchronization patterns which are asymptotically 
either optimal or almost optimal. In Section [Vl] we show 
how folding can be applied to construct multidimensional 
error-correcting codes. In section [VTIl we generalize the con- 
struction in [10] to form pseudo-random arrays on different 
multidimensional shapes. We conclude in Section rVlIII 

II. Folding and Lattice Tiling 

A. Folding 

Folding a rope, a ruler, or any other feasible object 
is a common action in every day life. Folding an one- 
dimensional sequence into a ZJ-dimensional array is very 
similar, but there are a few variants. First, we will summarize 
three variants for folding of an one-dimensional sequence 
sqSi ■ ■ ■ s,n-i into a two-dimensional array A. The general- 
ization for a /^-dimensional array is straightforward while 
the description becomes more clumsy. 
Fl. A is considered as a cyclic array horizontally and 
vertically in such a way that a walk diagonally visits 
all the entries of the array. The elements of the 
sequence are written along the diagonal of the r x t 
array A. This folding works if and only if r and t are 
relatively primes. 
F2. The elements of the sequence are written row by row 

(or column by column) in A. 
F3. The elements of the sequence are written diagonal by 
diagonal in A. 
Fl and F2 were used by MacWilliams and Sloane [10] to 
form a pseudo-random arrays. F2 was used by Robinson [8] 
to fold a one-dimensional ruler into a two-dimensional 
Golomb rectangle. The generalization to higher dimensions 
is straight forward. F3 was used in [9] to obtain a synchro- 
nization patterns in the square grid. 

B. Tiling 

Tiling is one of the most basic concepts in combina- 
torics. We say that a D-dimensional shape S tiles the D- 
dimensional space if disjoint copies of S cover . This 
cover of Z^ with disjoint copies of S is called tiling of Z^ 
with S. For each shape S we distinguish one of the points 
of S to be the center of S. Each copy of iS in a tiling has 
the center in the same related point. The set T of centers 



in a tiling defines the tiling, and hence the tiling is denoted 
by the pair {T,S). Given a tiling {T,S) and a grid point 
{ii, i2, ■ ■ ■ jio) we denote by c(ii, 12, ■ ■ ■ ,10) the center of 
the copy of S for which (ii, 12, ■ ■ ■ ,10) ^ S. We will also 
assume that the origin is a center of a copy of S. 

Lemma 1: For a given tiling {T,S) and a point 
(ii,i2, ■■■,iD) the point (ii,«2, ...,«!))- c(ii,i2, ••■,«£>) 
belongs to the shape S whose center is in the origin. 

One of the most common types of tiling is a lattice tiling. 
A lattice A is a discrete, additive subgroup of the real D- 



space 



W.l.o.g., we can assume that 



A = {-uiwi + U2W2 H yu^VD ■ ui, . . . e Z} (1) 



where {vi, W2, ■ ■ • , i^d} is a set of linearly independent vec- 
tors in R^. A lattice A defined by ([T]) is a sublattice of 
if and only if {vi, V2, ■ ■ ■ , vd} C 1P . We will be interested 
solely in sublattices of TP . The vectors wi, W2, . . . , are 
called basis for A C TLP , and the D y. D matrix 



G 



V21 



V12 
V22 



VDl VD2 



VlD 
V2D 

VDD 



having these vectors as its rows is said to be the generator 
matrix for A. 

The volume of a lattice A, denoted V'(A), is inversely 
proportional to the number of lattice points per unit volume. 
More precisely, V{K) may be defined as the volume of the 
fundamental parallelogram n(A) in MP , which is given by 



n(A)1i'{ei«i+6t'2- 



<DVD ■■Q<i^<l, ,l<i<D} 



There is a simple expression for the volume of A, namely, 
V{K) = |detG|. 

We say that A induces a lattice tiling of S if the lattice 
points can be taken as the set T to form a tiling (T, 5). 

III. The Generalized Folding Method 

In this section we will generalize the definition of folding. 
All the previous three definitions are special cases of the new 
definition. The new definition involves a lattice tiling (T, S), 
where S is the shape on which the folding is performed. 

A ternary vector of length D, {di, ^2, . . . , djj), is a word 
of length D, where di € {-1, 0, +1}. 

Let 5 be a Z?-dimensional shape and let S = 
{di,d2, ■ ■ ■ , do) be a nonzero ternary vector of length D. 
Let (T, S) be a lattice tiling defined by a ZJ-dimensional 
lattice A, and let S be the copy of S in (T, S) which includes 
the origin. We define recursively a. folded-row starting in the 
origin. If the point (ii , Z2, . . . , id) is in 5 then the next point 
on its folded-row is defined as follows: 

• If the point (ii + di, 12 + ^2, ••■,*£) + do) is in S then 
it is the next point on the folded-row. 

• If the point (zi + c?i, Z2 + ^2, ■ ■ • , *d + do) is in S' ^ 
S whose center is in the point (ci, C2, . . . , c^) then 

(ii + di- ci,i2 + d2 - 02, ■■■ ,iD + do - cd) is the 
next point on the folded-row. 



The new definition of folding is based on a lattice A, a 
shape S, and a direction S. The triple (A, S, S) defines a 
folding if the definition yields a folded-row which includes 
all the elements of S. It appears that only A and S determines 
whether the triple (A, S, S) defines a folding. The role of S 
is only in the order of the elements in the folded-row; and 
of course A must define a lattice tiling for S. 

How many different folded-rows do we have? In other 
words, how many different folding operations can be de- 
fined? There are 3^ — 1 non-zero ternary vectors. If A with 
the ternary vector {di,d2, . . . jdu) define a folding then also 
A with the vector {—di,—d2,...,—dD) define a folding. 
The two folded-rows are in reverse order, and hence they 
will be considered to be equal. Other than these pairs of 
folded-rows, we don't know whether for each D, there exists 
a f-dimensional shape S with ^ different folded-rows. 
An example for 13 = 2 is given next. 

Example 1: Let A be the lattice whose generator matrix 
given by the matrix 



G = 



3 2 
7 1 



One can verify that shapes tiled by this lattice have different 
folded-rows. It can be proved that this is the lattice with the 
smallest volume which has this property. 

The first two lemmas are an immediate consequence of 
the definitions and provide us a concise condition whether 
the triple (A, S, 6) defines a folding. 

Lemma 2: Let (T, S) be a lattice tihng defined by the 
D-dimensional lattice A and let 5 = {di,d2, ■ ■ ■ jdo) be 
a nonzero ternary vector. (A, S, S) defines a folding if and 
only if the set {{i-di,i-d2, ■■■ ,i-dD) — c{i-di, i-d2, ■■■ ,i- 
do) ■ < i < \S\} contains \S\ distinct elements. 

Lemma 3: Let {T,S) be a lattice tiling defined by the 
_D-dimensional lattice A and let 6 = (di, c?2, . . . , dri) be a 
nonzero ternary vector. (A, S, S) defines a folding if and only 
if (|5|-di, . . . , \S\-dD)-ci\S\-di, \S\-dD) = (0, . . . , 0) 
and for each i, < i < \S\ we have {i ■ di, . . . ,i ■ do) — 
c(i • di, . . . ,i • do) 7^ (0, . . . ,0). 

Before considering the general D-dimensional case we 
want to give a simple condition to check whether the triple 
(A, 5, 5) defines a folding in the two-dimensional case. For 
each one of the four possible folding definitions we will give 
a necessary and sufficient condition that the triple (A, 5, 5) 
defines a folding. 

Theorem 1: Let A be a lattice whose generator matrix is 
given by 

'" ^'ii V12 
V21 V22 



G = 



If A defines a lattice tiling for the shape S then the triple 
(A, 5, i5) defines a folding 

• with the ternary vector 5 ~ (+1,+1) if and only if 
g.c.d.(u22 - W21, wii - V12) = 1; 

• with the ternary vector 5 = (+1,-1) if and only if 
g.c.d.(u22 + W21, wii + V12) = 1; 

• with the ternary vector S = (+1,0) if and only if 

g.C.d.(wi2,W22) = 1; 



• with the ternary vector S = (0, +1) if and only if 
g.c.d.(t;ii,i;2i) = 1. 

Theorem [T] is generahzed for the Z3-dimensional case as 
follows. Let A be a D-dimensional lattice tiling for the shape 
S with the following generator matrix. 



Vjr + Vji for ^1 + 1 < r < f 1 + ^2, 



G = 



Vii Vi2 
V21 V22 



VlD 
V2D 



VDl VD2 ■ ■ ■ VdD 

Assume we have the direction vector 6 = 
(di, (i2, . . . , dzj). W.l.o.g. we assume that the first £i > 1 
values are H-l's, the next £2 values are -I's, and the last 
D — £1 — £2 values are O's. There exist integer coefficients 
ai, a2, . . . , q;d such that 

D 

^aj{vji,vj2,...,VjD) 

= (|5|,...,|5|,-|5|,...,-|5|,0,...,0), 

and there is no integer i,0 < i < \S\, and integer coefficients 

f3i,l32, . . . ,/3_D such that 



D 



Hence we have the following D equations: 

D 

ajVjr — \S\, I < r < £1, 



-1,0,. ..,0) 



D 



^ajVjr^-\S\, £i + l<r<£i+£2, 

D 

ajVjr = 0, £i+£2 + l<r < D, 
which are equivalent to the following D equations: 

D 

^ttj-Wji = \S\, 



D 



ctjivjr — Vji) = 0, 2 < r < ^1, 



aj{vjr + Vji) = 0, £i + l <r <£i+£2, 



D 



Y "i^jr = 0, £i+£2 + l<r < D. 
i=i 

We define now a set of _D(Z? — 1) new coefficients Urj, 
2 <r < D, 1 < j < D, CIS follows: 



for £1 



Kr < D. 



Consider the {D — 1) x D matrix 



H 



U21 U22 
U31 U32 

UDl UD2 



U2D 
U3D 

udd 



v,i for 2 < r < £1, 



Using the Cramer's rule it is easy to verify that the 
following assignment is the unique solution for the ai's. 

= (-l)'"MetiJ„ l<i<D, 

where Hi is the (Z) — 1) x (Z? — 1) matrix obtained from H 
by deleting column i of H. 

Theorem 2: If A define a lattice tiling for the shape S 
then the triple (A, 5, 5) defines a folding if and only if 
g.c.d.(dctiJi,dctiZ2, ■ • . ,deti7i3) = 1. 

One important tool that we will use to find an appropriate 
folding for a shape S' is to use a folding of a simpler shape 

5 with the same volume and apply iteratively the following 
theorem. 

Theorem 3: Let A be a lattice, S a D-dimensional shape, 

6 = (di, ^2, . . . , dn) a direction, and {A,S,S) defines a 
folding. Assume the origin is a point in S, (ii, «2, • • ■ , *l>) G 
S, (ii + di,i2 + d2, . . . ,iD + do) £ S, S ^ S, and the 
center of S is the point {ci,C2, ■ ■ ■ ,cd)- Then for the shape 
S' = S U {(ii + di,i2 + d2, ■ ■ ■ ,iD + do)} \ {{h + di - 
ci,i2 + d2-C2, ■ ■ ■ ,iD+dD~CD)} the tiiple (A, 5', (5) also 
defines a folding. 

IV. Bounds on Synchronization Patterns 

Our motivation for the generalization of the folding op- 
eration came from the design of two dimensional synchro- 
nization patterns. Given a grid (square or hexagonal) and 
a shape S on the grid, we would like to find what is the 
largest set V of dots on grid points, \V\ — m, located in 
S, such that the following property holds. All the (™) lines 
between dots are distinct either in length or in slope. Such a 
shape S with dots is called a distinct difference configuration 
(DDC). If 5 is an TO X TO array with a dot in each row and 
a dot in each column than S is called a Costas array [5]. 
If 5 is a fc X TO array with a dot in each column then S 
is called a sonar sequence [5]. If iS is a fc x n array then 
S is called a Golomb rectangle [7]. These patterns have 
various applications as described in [5]. A new application of 
these patterns to the design of key predistribution scheme for 
wireless sensor networks was described lately in [12]. In this 
application the shape S might be a Lee sphere, an hexagon, 
or a circle, and sometimes another regular polygon. This 
application requires in some cases to consider these shapes 
in the hexagonal grid. F3 was used for this application in [9] 
to form a DDC whose shape is a rectangle rotated in 45 
degrees in the square grid. Henceforth, we assume that our 
grid is the square grid, unless stated otherwise. 

Let S and S' be two-dimensional shapes in the grid. We 
will denote by A{S,S') the largest intersection between S 



and S'. \S\ will denote the number of grid points in S. Let 
TO be a given integer An infinite set of dots in the grid such 
that each given shape S on the grid is a DDC with to dots 
will be called an infinite S-DDC. The following theorems 
are generalization of similar theorem in [9]. 

Theorem 4: Assume we are given an infinite 5-DDC with 
TO dots on the grid. Let TZ be another shape on the grid. Then 
there exists a copy of TZ on the grid with at least A(5, TV) 
dots. 

Theorem 5: Assume we are given an infinite 5-DDC with 
TO dots on the grid. Let Tl and U be another shapes on the 
grid. Then there exists a copy of U on the grid with at least 
l;5^A(5,7^) • l\{n,U) dots. 

In order to apply Theorem |4] and Theorem |5] we will use 
folding of the sequences, defined as follows, in our shape S. 
Let A be an abelian group, and let E = {ai, 02, ... , am] C 
A be a sequence of to distinct elements of A. We say that 
i? is a B2-sequence over A if all the sums a,;j + a^^ with 
1 < ii < i2 < ni are distinct. For a survey on _B2-sequences 
and their generalizations the reader is referred to [13]. The 
following lemma is well known and can be readily verified. 

Lemma 4: A subset E = {ai, a2, . . . , C A is a B2- 
sequence over A if and only if all the differences a^^ — a^^ 
with 1 < «! 7^ 12 < ™ are distinct in A. 

Note that if D is a i32-sequence over Z„ and a G Z„, 
then so is the shift a + E = {a + e:e^ E}. The following 
theorem, due to Bose [14], shows that large i32-sequences 
over Z„ exist for many values of n. 

Theorem 6: Let g be a prime power Then there exists a 
_B2-sequence oi, 02, . . . , a,„ over Z„ where n = — 1 and 
m = q. 

The importance of folding a B2 sequence S into a given 
shape S is given by the following theorem. 

Theorem 7: Let A be a lattice, S, n — \S\, a D- 
dimensional shape, and S a direction. Let E be a B2- 
sequence over Z„. If {A,S,S) defines a folding then the 
folded-row is a /^-dimensional DDC. Moreover, this DDC 
can be extended to infinite 5-DDC. 

In the sequel we will use Theorem |4] Theorem |5] and 
Theorem |7] to form DDCs with various given shapes with a 
large number of dots. To examine how good are our bounds 
on the number of dots we should know what is the upper 
bound on the number of dots in a DDC whose shape is S. It 
was shown in [9] that for a DDC whose shape is a regular 
polygon or a circle, the maximal number of dots is at most 
^/i + o(^/s), when the shape contains s points of the grid. 

V. Bounds for Specific Shapes 

In this section we will present some lower bounds on the 
number of dots in some two-dimensional DDCs with specific 
shapes. One of the main keys of our constructions, and the 
use of the given theory, is the ability to produce a DDC with 
a rectangle shape and any given ratio between its edges. 

Theorem 8: For each positive number 7 there exist two 
integers a and b such that | w 7 and an infinite 5-DDC 
with Va • bR + o{R) dots whose shape is an a x /3 = {bR + 
o{R)) X (aR + o{R)) rectangle, a(3 = — 1 for some prime 
p, and even a. 



Now, we can give a few examples for specific shapes. To 
have some comparison, let the radius of the circle or the 
regular polygons be R (the radius is the distance from the 
center of the regular polygon to any one its vertices). 

A. Regular Hexagon in the Square Grid 

By Theorem [8] there exists an infinite 5-DCC, where S 
is an a X /3 = [aR + o(i?)) x {bR + o{R)) rectangle, such 
that ^ ~ ctP = p^ — 1 for some prime p, and even a. 
Let A be the a lattice tiling for S with the generator matrix 

a 

where 9 — 1 if a = {mod 4) and 6 = 2 if a = 
2 (mod 4). By Theorem [T] (A, 5, J), S = (+1,0), defines 
a folding. Now, Theorem [3] is used iteratively to form an 
infinite 5'-DCC, where S' is "almost" a regular hexagon 
with \/ a ■ bR + o{R) dots (the six vertices of one hexagon 
are at (f , 0), 0), (f , f ), a), (f , a), (0, f )). Hence, 
a lower bound on the number of dots in a regular hexagon 

with radius R is approximately R + o{R). 

B. Circle in the Square Grid 

We apply Theorem |4] where 5 is a regular hexagon with 
radius p and S' is a circle with radius R. A low er bound on 

the number of dots in S is approximately p+o{p). The 

maximum on ^^^^^^° —A{S, S') yields a lower bound of 
1.70813R + o{R) on the number of dots in S'. 

C. Other Shapes 

For most of the regular rt-gons {n ^ {4, 6}) in the square 
grid we applied Theorem |5] with a hexagon and a circle 
(for n ~ 4 the optimum can be obtained from a Costas 
array). Exceptions are n = 3, 5, 8, 10, and 12, where 
we got a better bound by specific constructions. Table H] 
summarize the bounds we obtained for regular polygons and 
a circle in the square grid. We also consider circle in the 
hexagonal grid, but the main result in the hexagonal grid is a 
construction of an optimal DDC whose shape is a hexagon. 
This involves another interesting construction for a lattice 
tiling of another shape in the square grid and translation into 
the hexagonal grid. The same techniques can be used for any 
_D-dimensional shape. Finally, we note that the problem is 
of interest also from discrete geometry point of view. Some 
similar questions can be found in [15]. 

VI. Application in Error-Correction 

Assume that we have a D-dimensional array of size ni x 
7i2 X • • • X and we wish to correct any D-dimensional 
burst of length 2 (at most two adjacent positions are in error). 
The following construction given in [16] is based on folding 
the elements of a Galois field with characteristic 2 in a parity 
check matrix, where the order of the elements of the field 
is determined by a primitive element of the field. 
Construction A: Let a be a primitive element in GF(2™) 
for 2™ - 1 > Hfci '"'i- Let d = {log^ D'] and i = 
(ii, Z2, . . . , in), where < < — 1. Let Ahe a d x D 



TABLE I 

Bounds on the number of dots in an n-coN DDC 



n 


upper bound 


lower bound 


ratio between bounds 


3 


1.13975-R 


1.02462_R 


0.899 


4 


1.41421_R 


1.41421_R 


1 


5 


1.54196-R 


1.45992_R 


0.9468 


6 


1.61185-R 


1.61185-R 


1 


7 




1 p;c;9'5QR 

1 .OOZOO-fl. 


u.yoo^ 


8 


1.68179R 


1.60094_R 


0.9519 


9 


1.70075ii 


1.61343iJ 


0.9487 


10 


1.71433ii 


1.64302_R 


0.9584 


11 


1.72439i? 


1.64458i? 


0.9537 


12 


1.73205iJ 


1.66871i? 


0.9634 


13 


1.73802_R 


1.66257iJ 


0.9566 


14 


1.74275i? 


1.66883iJ 


0.9576 


30 


1.76598iJ 


1.69955iJ 


0.9593 


90 


1.77173JJ 


1.70718iJ 


0.9636 


circle 


1.77245R 


1.70813R 


0.9637 



matrix containing distinct binary c?-tuples as columns. We 
construct the following ni x n2 x • • • x no x (m + d + 1) 
parity check matrix H. 



hi = 



1 

mod 2 



for all i = («!, «2, ■ • ■ , io), where < ig < ni — 1. 

Theorem 9: The code constructed in Construction A can 
correct any 2-burst in an rti x rt2 x • • • x array codeword. 

Theorem 10: The code constructed by Construction A has 
redundancy which is greater by at most one from the trivial 
lower bound on the redundancy. 

The same construction will work if instead of a D- 
dimensional array our codewords will have have a shape S 
of size 2™ — 1 and there is a lattice tiling A and a direction 
5 such that (A, 5, 5) defines a folding. The elements of 
GF(2™) will be ordered along the folded-row of S. The 
construction can be generalized for more complicated types 
of multidimensional errors. 

VII. Application in Pseudo-Random Arrays 
Let n = 2^^^'' - 1 such that ni = 2^^ - 1 and n2 = — ai-e 

^ ^ ni 

relatively primes and greater than 1. Let S = sqSi ■ ■ ■ s„_i 
be an m-sequence (maximal length linear shift register 
sequence [10], [17]) of length n. Assume we use Fl to form 
an m x n2 array A. A has many interesting properties such 
as shift, recurrences, addition, auto-correlation, etc. [10]. It 
also has a fci x fc2 window property, i.e., each fci x ^2 
possible binary matrix appears exactly once as a window 
in the cyclic array. These arrays were called in [10] pseudo- 
random arrays. All these properties except for the window 
property are a consequence of the fact that the elements in 
the folded-row are consecutive elements of the m-sequence 
S. Hence, if we use any of the folding operations to fold 
S into a D-dimensional shape S, the shape S will have all 
these properties. As a consequence of Theorem [3] we have 
the following theorem. 



Theorem 11: Assume A define a lattice tiling for an ni x 
7i2 array, such that ?iin2 = 2^^^'^ — 1, ni, rt2 are relatively 
primes and greater than 1. Assume further that A defines a 
lattice tiling for the shape S and (A, 5, 5) defines a folding 
for a direction 5. Then, if we fold an m-sequence S into S 
by the direction 5 then the resulting shape S has the ki x k2 
window property if and only if the ni x n2 array A has 
the ki X fc2 window property by folding S into A by the 
direction 5. 

VIII. CONCLUSION 

The well-known definition of folding was generaUzed. 
The generalization makes use of a lattice tiling and a di- 
rection in which the folding is performed. We demonstrated 
how folding in general and the new definition in particular 
is applied for constructions of multidimensional synchro- 
nization patterns, error-correcting codes, and pseudo-random 
arrays. The compressed discussion we have made raised 
lot of problems for further research, which can advance 
the research on multidimensional coding. It also raised 
intriguing questions which are related to discrete geometry. 
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